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I. INTRODUCTION 

This Final Report covers the work done by ACM under 
Contract NAS9-15445 during the time period November 1, 1977 
to March 1, 1979. This report will serve as an overview and 
summary of the work done. More details on the work are avail- 
able from the individual reports listed in the references. 

1.0 Orbit Computation \u>J hods 

The methods for comput ing satellite* orbits can generally 
be put into two classifications: 

(1) Nurr,.: fi e.j ; Uc - The solution to the satellite 

differential equations are solved in a step-by-step 
manner using a mathematical algorithm taken from 
numerical analysis. Very precise orbits can be 
computed in this way. Detai led models o! the per- 
turbing forces are allowed. An important disadvan- 
tage is that these methods usually require signifi- 
cant amounts ol computer execution time. Hence 
they are impractical tor many apolications where 
rapid calculations are required. These methods are 
also referred to as cp. ’ial perturbations . 

(2) nnalytical Methods - The solution is expressed by 
explicit functions of the independent variable. 
Evaluation ol these functions yield the state vec- 
tor of the satellite at the desired time. No inter- 
media 1 •• calculations are necessary. Accordingly, 
rapid generation of satellite enhemerides is possi- 
ble. These methods are also referred to as general 
pertu rba t ions . 
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Tilt' various available methods have boon applied to astro- 
dynamics problems, according to their particular advantages. 
However, analytical methods have not found wide acceptance 
because they usually did not contain sufficiently precise 
models of the per turbine forces. Also, a large initial invest 
nient in manpower is required in order to develop the explicit 
functions. These disadvantages have been overcome in recent 
years due to advances in celestial mechanics. The purpose of 
the work described in this report was to apply an up-to-date 
analytical orbit calculation method in the development of an 
efficient and accurate orbit prediction program. 


2 . 0 App 1 i ca 1 i ons o f A na 1 y t i ca 1 '!< »thcn i_e 

Wit), the development of NASA's Space Transportation Sy- 
stem, ana in particular the Shuttle Orbiter, new demands are 
being made on computer software for orbit computation. The 
emphasis is on rapid calculations in order to support the 
many planned STS flights. Also, t ho increased complexity of 
STS flights means that the overall computation load will be 
increased, requiring that the most efficient mathematical 
algorithms be used in t he software. For t hose reasons, 
analyt ical methods arc preferred for many applicat ions in 
support of Shuttle missions. They are particularly appli- 
cable to iterative calculations such as targeting routine's 
and orbit determination routines, and in other cases where 
time critical orbit calculations are required. 
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1.0 Dackgreui i d 

A fundamental new idea in celestial mechanics was pre- 
sented by Dr. Gerhard Scheil'ele in 1970. lie expressed the 
satellite equations of motion in an extended phase space, 
us ini; the true or eccentric anomaly as the independent vari- 
able. Since the formulation was canonical, it touiul direct 
application to analytical satellite theories. 

A set o I orbital e 1 * mu Mils similar t « > Delaunay's < • I * mien Is 
(cal led DS-e I omen t.s ) were developed by mu kill!'; us** ol tin* new 
I oriiia 1 t sin . It was round that those elements presented signif- 
icant advantages in developing analytical theories. A solu- 
tion to the J., oblateness problem was carried out, and it was 
shown that this solution was more accurate than solutions 
based on classical theories. 

Singularities at zero eccentricity and inclination were 
eliminated by the development of yet another set ol extended-space 
canonical elements: the Poincare-Similar elements in true 

anomaly, called PSg-e lemon ts . The J,, solution in true anomaly 
DS-elements was converted to PSi*>-elements . 

Major advances to the realization of a completely ana- 
lytical solution to the near-earth satellite problem were 
made by ACM under NASA/JSC Contract NAS9-15171. For the first 
time, realistic models of the perturbing forces were included 
in the solution. The solution o! a drag-perturbed satellite 
was carried out, using an atmospheric density model that in- 
cludes all the important periodic variations in density. In 
addition, the implementation of this model into the' analytical 
solution allows for coupling between drag and J., , and between 
drag and itself. The long period effects due to the zonal 
geopotential terms were also added to the solution. An arbi- 
trary number of zonal terms is allowed. 



Also, under Contract NAS9-15171, the first steps were 
made toward tin* conversion ol the analytical solution into 
etticient computer software. A modular structure was 
developed in order to allow maximum flexibility in applica- 
tions. The short period terms of the solution were induced 
in the software, that was called the Analytical Satellite 
Orbit Prediction Program (ASOP). 


2.0 Scope 

Under Contract NAS9-15445, further developments were 
carried out on the ASU'* computer program. (The previous 
version of ASOP contained only the short period J,, terms.) 
Particular emphasis was given to building a modular program 
structure. The work can be divided into the following major 
parts : 

(1) Furt her Mathematical Development - Analytical 
initialization formulas were developed for includ- 
ing the effects of the tosseral (time dependent) 
geopotential terms. Further refinements were 
made on the drag solution, with the goal of estab- 
lishing a concise formulation. 

(2) Development of Program Structure - The modular pro- 
gram structure was developed, including module 
interlaces and data management. 

(3) Development of Program Modules - Modules were 
developed for (a) atmospheric drag, (b) tesseral 
geopotentin 1 terms, (c) «!„ short period and first 
order secular terms, (d) J 2 second order and long 
period terms, and (e) second order secular and long 
period terms due to higher order zonals. The com- 
puter code of each module was optimized with respect 
to computer storage. 

(4) Program Documentation - Program documentation and 
user's guide was prepared for ASOP. 
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(5) Pro^ran Pert wd Verification - Individual modules 
and options w< • tested. The combined program was 
verified by comparing to a precision numerical inte- 
gration. 

All modules that were developed in this work make use of the 
PSO- elements that were developed by ACM under previous NASA 
contracts. These are Poincare-Similar elements based on the 
true anomaly as independent variable. Therefore, the ASOP 
propram contains no singularities for zero • •centricity or 
inclination. The modular structure o f ASOP allows the user 
to build an analytical orbit prediction program that is con- 
sistent with his computer memory restrictions and accuracy 
requirements. The ASOP program documentation is given in 
Reference 1. 
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III. SUMMARY OF THE INDIVIDUAL TASKS 


1.0 Analytic Prut; Module 

A set of subroutines was developed that gives the atmos- 
pheric drag effects on a satellite orbit. They were implemented 
into ASOP. 

The atmospheric density model and drag solution were taken 
from the analytical developments carried out under Contract 
NAS9-15171. The module code was compacted and optimized with 
respect to computer run* ime and storage*. More details on the 
work are given in References 1 and 2. 


2.0 Tessera 1 Terms Initialization Module 

The dominant perturbations of the motion of a satellite 
near the earth are due to the non-symmetrical gravitational 
field and the atmospheric drag. The gravitational field may 
be divided in two classes: terms independent of time (zonal 

harmonics) and terms which depend explicitly on time (tesseral 
harmonics). This task concerns the unified treatment of the 
tesseral terms in the PS$ theory and their implementation 
into ASOP. 

The perturbations due to the tesseral harmonics can be 
placed in four categories: 

(1) Short period perturbations with a magnitude of 
about J 2 

(2) Intermediate period perturbations with a magnitude 
of between J n and J, 2 , 

(3) Along track secular perturbations induced by the 
periodic perturbations in the mean motion. 

(4) Resonant perturbations. 

For near earth satell ites , the atmospheric drag perturba- 
tion continually pulls the satellite in and out of the different 
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long period resonant frequencies. The result is that the 
resonances never become app rent and may be neglected. 

Since the J 9 theory has been deve loped only to first, 
order, the second order short period tessera 1 perturbations 
may be neglected. For the same reason, the intermediate pe- 
riod perturbations should bo included if first order accuracy 
is to be maintained. The tessera 1 hnrn.onic.s have no true sec- 
ular perturbation but the poriodieit ies in the 1 mean motion 
induce a secular perturbation in the mean anomaly. This sec- 
ular perturbation may be determined by simply using the aver- 

a. 

age mean motion instead of the osculating mean motion. Graf 
finds the average moan motion in a numerical manner. The nu- 
merical studies with that approach verify the assumption that 
use of the average mean motion does account for the apparent 
secular trend in the mean anomaly. Although the results are 
good, use of a numerical method would be inconsistent with 
the idea of having a completely analytical theory. 

To complete the solution of the motion of a near earth 
satellite, the averaged mean motion and the intermediate pe- 
riod perturbations wore found in a completely analytical 
manner . 

Since the previous l\S if) theory applied Von Zeipel's solu- 
tion technique, it seemed natural to return to this method 
for a solution of the tesseral perturbation. As in the pre- 
vious developments, the solution was first, developed in DS'|j- 
elements, then converted to PS<{)- elements to avoid singular- 
ities. Details of the work are given in Reference 3. Numer- 
ical testing of this module is described in Reference 4. 


See Final Report on Contract NAS!) I . r > I 7 1 . 
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3.0 Second Order Secular and I >ng Period Terms Modulo 

The second order secular and long period terms due to 
all zonal harmonies were implemented into ASOP. The analyt- 
ical formula i ion was taken from the work don*' under Contract 
NAS9-15171 . 

This module contains two parts: 

(1) Terms due to and, 

(2) Terms due to all necessary higher order zonal 
terms . 

(The work done in building this module covered two tasks in 
Contract NAS9- 15445 Statement of Work). Documentation of the 
equations in this module are given in Reference 5. A report 
on the testout of this module is given in Reference 6. 


4.0 Verification Testing of t he A SOP Program 

The accuracy of each analytical theory implemented in 
the ASOP program was determined during the development of the 
model. however , a final eheekoul was requ i red i<> further de- 
termine the accuracy and the limits of the whole program after 
all the models had been included. Therefore, a series of 
tests were designed in order to evaluate the total program. 
Twenty- five (25) test cases were run on a UNI VAC 1110 system. 

A detai led descripi ion of tesl cases and discussion of re- 
sults are given in Reference 7. The analytical solutions 
obtained f rom ASOP wore compared to numerical solutions 
obtained from the KSFAST program. 
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